Scars on quantum networks ignore the Lyapunov exponent 
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We show that enhanced wavefunction localization due to the presence of short unstable orbits and 
strong scarring can rely on completely different mechanisms. Specifically we find that in quantum 
networks the shortest and most stable orbits do not support visible scars, although they are respon- 
sible for enhanced localization in the majority of the eigenstates. Scarring orbits are selected by 
a criterion which does not involve the classical Lyapunov exponent. We obtain predictions for the 
energies of visible scars and the distributions of scarring strengths and inverse participation ratios. 
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One of the most striking ways in which the underly- 
ing classical dynamics of a chaotic system manifests it- 
self in the correspondingquantum behavior is the scar 
phenomenon 0, H H UM 0, II II H E3 A scar is a 
quantum eigenfunction with excess density near an un- 
stable classical periodic orbit (PO). Such states are not 

i expected within Random-Matrix Theory (RMT), which 
predicts that wavefunctions must be evenly distributed 
over phase space, up to quantum fluctuations Ex- 
perimental evidence and applications of scars come from 

[ systems as diverse as microwave resonators |^, quantum 
wells in a magnetic field 0, Faraday waves in confined 

[ geometries Id, o pen quantum dots '8!| and semiconductor 
diode lasers jinl |. 

Quantum networks (graphs) are established models in 
the field of mesoscopic physics, from which most of the 
above examples are drawn, as well as in many other areas 

' including molecular and mathematical physics and quan- 
tum computation (see 0,0,0,0 and Refs. therein). 
In recent years they have become one of the most promi- 
nent tools in quantum chaos because they allow to study 
with simple means the applicability of RMT and its limi- 
tations due to system-specific properties 0, 0, 0] . For 

I example, it was shown recently [sj that statistics of the 
bulk of graph eigenfunctions (including, e. g. the left 
eigenstate in Fig. ^) conform with the existing theories 
describing the effect of short unstable periodic orbits on 

' the localization properties of wavefunctions. 

Therefore it is surprising that the same does not ap- 
ply to the small but important group of strongly scarred 
eigenstates (Fig. right), which we study in this let- 
ter. We show that, contrary to common intuition and 
accepted theories 0, 12 S S 0|j the shortest and least 
unstable orbits of the system produce almost no visible 
scars, although they are responsible for enhanced local- 
ization within the bulk of states. We derive a condition, 
Eq. O below, selecting orbits relevant for strong scar- 
ring from topological information only and without any 
reference to the classical Lyapunov exponent. Based on 
this insight we are able to give a criterion for the energies 
at which strong scars are to be expected, Eq. iPJ, and de- 



scribe their statistical distribution, Eq. (|12|l . In view of 
the numerous and diverse applications of quantum net- 
works 0, 0, 0, 0, 0| , these main results should be 
of broad interest in their own right. On top of this, some 
important conclusions generalize beyond graphs. In par- 
ticular, our results provide clear evidence for the fact 
that enhanced wavefunction localization due to the pres- 
ence of short unstable orbits and strong scarring can in 
principle rely on completely unrelated mechanisms and 
can also leave distinct traces in statistical measures such 
as the distribution of inverse participation ratios (IPR). 

Following the quantization outlined in |l4l | we consider 
graphs which consist of V vertices connected by a set of B 
bonds. The number of bonds emanating from a vertex j 
is its valency Vj . A basis state on the network is specified 
by a directed bond d = [i ^ j], i. e. an ordered pair of 
connected vertices i, j. Hence a quantum wavefunction 
is just a set of 2B complex amplitudes Ud, normalized 
according to J2d\^d\'^ ^ ^- "^^e standard localization 
measure is the IPR 

2B 

i=Y.\^d\'. (1) 

d=l 

Ergodic states which occupy each directed bond with the 
same probability have 2 = 1/2B and up to a constant 
factor depending on the presence of symmetries this is 
also the RMT prediction. In the other extreme T = 1 
indicates a state which is restricted to a single bond only, 
i. e. the greatest possible degree of localization. 

Fig.n^ shows the distribution of J for fully connected 
graphs. Some features of this distribution are explained 
by the original scar theory of Heller and extensions of 
it 0, 1^ . The main idea is to connect localization proper- 
ties of eigenfunctions to the dynamics of the system. For 
example, the identity {T) = limT^oo /o'"dt(P(t)) ex- 
presses the mean IPR in terms of the quantum return 
probability (RP) P{t), averaged over time and initial 
states. It is then argued that the short-time dynamics, 
approximated semiclassically with some short PO's, pro- 
vides sufficient information for estimating (Z). Within 
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FIG. 1: Statistical properties of eigenstates for fully con- 
nected graphs with different valencies v = V — 1. The eigen- 
states were obtained by diagonalization of the bond-scattering 
matrix Eq. Q, and a statistical ensemble was generated by 
choosing random bond lengths, (a) The bond intensities of 
two representative eigenstates are shown with a gray-scale for 
a typical state (left) and a scar on a triangular PO (right), (b) 
Probability distribution of the IPR, showing a step-like cutoff 
at X = 1/6 that can be attributed to scarring on triangular 
orbits, (c) The bulk of the IPR distribution shows scaling ac- 
cording to Eq. ((21 . (d) The quantum return probability and a 
classical approximation (*) based only on the period- two or- 
bits; (e) Rescaled and integrated distribution of triangle scars 
Vint{S) = B-^ dS'P^^\S'). The dashed line has slope 1, 
corresponding to the theoretical prediction Eq. 1)12^ . 



this approach it is clear that the orbits with the low- 
est Lyapunov exponent (LE) have the largest influence 
on eigenfunction localization because classical trajecto- 
ries can cycle in their vicinity for a relatively long time 
and increase the RP beyond the ergodic average. For the 
graphs studied in this letter, some PO's p and their LE 
Ap are listed in Table |l| The shortest PO's have period 
2 and bounce back and forth between two vertices. For 
large graphs i; ^ oo these are by far the least unstable 
ones, as their LE approaches while all others become in- 
creasingly unstable Ap ~ Inw. Indeed the period 2 orbits 
totally dominate the classical and quantum RP at short 
times (see Fig. QJi). Including the contribution of these 
orbits only, Kaplan obtained a mean IPR which is by a 
factor V larger than the RMT expectation, in agree- 
ment with numerics 0|. Moreover, following the same 
line of argumentation as in Q we get that the bulk of 
the IPR distribution scales as 



TABLE I: The topology of the shortest PO's of a fully con- 
nected graph with valency v are shown with the corresponding 
amplitude, Lyapunov exponent, and IPR. 
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indicating that the whole bulk of P(I) is effectively de- 
termined by the least unstable orbits (Fig.^). 

With all this evidence for their prominent role in wave- 
function localization, one clearly expects to see strong 
scarring on the period 2 orbits [IJI ■ Such states would es- 
sentially be concentrated on two directed bonds and give 
rise to X ~ 1/2. However, in this region V{I) is negligi- 
ble. The shortest and least unstable orbits of our system 
produce no visible scars. Note that the same applies also 
to the value X = 1/4 expected from the V-shaped orbits 
of Tabled In fact V{I) has an appreciable value only 
for Z < 1/6 (Fig. nja). The position of this cutoff pre- 
cisely coincides with the IPR expected for states which 
are scarred by triangular orbits. They occupy six directed 
bonds since, due to time-reversal symmetry, scarring on 
a PO and its reversed must coincide. Indeed a closer in- 
spection shows that the vast majority of states at Z « 1/6 
look like the example shown in Fig.Q^ (right). Of course 
the step at Z = 1/6, which is present for any graph size V, 
is incompatible with the scaling of P{I) mentioned above 
and indeed this relation breaks down in the tails at the 
expected points (inset of Fig.QJ;). We conclude that visi- 
ble scars on short unstable orbits can strongly modify the 
tails of the IPR distribution even beyond the known pre- 
dictions for the influence of short PO's on wavefunction 
localization. 

In the rest of the paper we will formulate a theory 
which explains the above observations. To this end we 
must be more explicit concerning the dynamics on the 
graph. We consider particles with fixed wavenumber k, 
propagating on the bonds and scattering at the vertices. 
During the free propagation on a directed bond the wave- 
function accumulates a phase fciij, where Lij = Lji 
denotes the length of the corresponding bond. At the 
vertices current conservation and continuity of the wave 
function are required. These boundary conditions can be 
translated into vertex scattering matrices, which describe 
a unitary transformation from Vi incoming to Vi outgoing 
waves at each vertex. Without loss of generality we re- 
strict the presentation to the simplest case of Neumann 
boundary conditions jl^ , where the scattering matrix of 
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vertex i is 



2/v, 



(3) 



reasoning that lead to Eq. © now to the bonds D^- J and 
making use of Eq. Q we get 



We can now combine the free propagation and the vertex 
scattering into a 2B x 2B operator, the bond-scattering 
matrix S 14] , which acts on the amphtudes associated 
with the directed bonds. The matrix element 



e'^ ad 



{d e D) 



(8) 
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which relates the amplitude on a directed bond d E p to 
the amplitude on the reversed bond d G p. This means 
that the same states are scarred on p and p, as expected 
from time-reversal symmetry. Substituting d ^ d we get 



describes a transition from the directed bond d — [i j] 
to d' = [m — > n]. We interpret S as quantum time- 
evolution operator on the graph. {S^)d'd is the complex 
probability amplitude to be after t — 0,1, . . . time steps 
on the directed bond d' if the initial state was on d. In 
particular, \{S'^)dd\'^ is the quantum RP shown in Fig.^. 

.ikL 



ikL 



Qd which together with Eq. JHJ implies 



(kLd — A) mod vr = y d £ p 



(9) 



with arbitrary A (18i] . Eq. (^J is a necessary and sufficient 
condition for the energies of perfect scars. It is reminis- 
cent of a simple Bohr-Sommerfeld quantization condition 



= Ep Ape' " is expanded as a sum over aU PO's ^ 2mT + A, as it applies, e. g., to strong scars in bil 



of period t starting at d. Here Lp is the total length of 
orbit p. Assuming for simplicity vj = v throughout the 
graph, we express the amplitude Ap by the number rp of 
reflections along p, Ap = (2/t;)*-''f (2/u-l)''p (cf Table|Ill . 
The classical RP is obtained by summing, instead of the 
amplitudes, the probabilities Mp = As Mp < 1, 

the probability to follow the PO decreases exponentially 
with time and the orbit is unstable. Hence one defines 
the Lyapunov exponent of a PO p on a graph by 



A^ 



-tpHnMp. 



(5) 



Let us now come back to the problem of scarring and 
investigate the conditions under which we can construct 
perfect scars on the graph, i. e. eigenstates S\a) = e'^|a) 
which have the property that they have non-zero am- 
plitude only on a PO p and vanish on all other bonds. 
Consider an arbitrary vertex j on p and let -D^ be the 
set of directed bonds which are leaving/ approaching j 
and belong to p. Similarly let d'-'^ be the set of bonds 
not belonging to p. By construction, the amplitude of \a) 
on these latter bonds vanishes, i. e. all waves arriving on 
the bonds d'j~J and transmitted across the vertex to a 
bond in D^^J cancel each other 



d' 



dd'Od' = — 
Vi 
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ikL 



^' ad' (den'^i^). (6) 



This equation has an important consequence: a perfect 
scar cannot live on a single bond attached to j because 
there would be no way to cancel the transmitted wave. 
An exception are only vertices with valency Vj = 1, for 
which d'^^ is empty. Formally a necessary condition for 
scarring orbits p is 



> 2 - 



(Vj e p) 



(7) 



where Vj^p denotes the number of bonds attached to j 
and belonging to p. Eq. Q excludes, in particular, per- 
fect scars on the period-two orbits. Applying the same 



liards. However, there is an important difference: not 
only does Eq. ^ require quantization of the total action 
kLp of the scarred orbit, it also implies action quantiza- 
tion on all the visited bonds d. This stronger condition 
can only be met if the lengths of all bonds on p are ratio- 
nally related. As in general the bond lengths are incom- 
mensurate there are no perfect scars for generic graphs. 

Nevertheless, for incommensurate bond lengths Eq. @ 
can be approximated with any given precision and then 
visible scars are expected. A natural measure for the 
quality of a scarred state \a) is the total probability Sp — 
'^d^p l^rfP to find this state away from the scarring orbit 
{Sp = corresponds to a perfect scar). We will derive the 
probability density of strong visible scars P{Sp — > 0). We 
represent the bond-scattering matrix S as perturbation 
of a matrix Sq which has a perfect scar on p, i.e. 



S^e+'"^ So « (l + ie$)S'o. 



(10) 



The deviations from exact quantization for the individ- 
ual bonds d € p have been combined into a diagonal 
matrix $ with e^dd — {kLd ^ A) modTr for d (z p,p and 
^dd = otherwise. The strength of the perturbation is 
given by Sp = minx max^gp \ {kLd — A) mod7r|. For a PO 
covering N undirected bonds of the graph N bond phases 
kLdiaodir must approximately coincide. Upon variation 
of k, they are independent and uniformly distributed ran- 
dom numbers in [— tt, +7r] . Therefore the probability den 
sity of a small perturbation is p{e — *■ 0) 
first order perturbation theory we have 



a'^^^\^^>\a'^'>){d\a)^>) 



(0)\ 



(11) 



where A 



(0) i.(o)\ 



are eigenphases and eigenvectors of 5*0 , 
including the perfect scar |al°''). The quantity x = Sje^ 
is distributed with some probability density p{x) that 
is independent on e. Consequently we have p{S\£) = 
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FIG. 2: ViT) for star-graphs of various valencies (inset a 
star-graph with i; = 6). In the upper inset we report the 

integrated distribution 'Pint (5) = /p dS'P{5') in a double- 
logarithmic plot. The dashed line has slope 0.5 which cor- 
responds to the theoretical prediction Eq. I12II . 

e~^p((5/e^). We can now use 'P{S) — J de p{S\e) p{e) to- 
gether with p{e ^ 0) e^^^ to deduce 

=CJ(^-3)/2 (^^0). (12) 

Here C = J dx x^'^^^^''^^ p{x) is a constant which de- 
pends on the size and topology of the graph. Note that, 
according to Eq. x is a sum oi B — N independent 
non-negative terms x = X^d^p \xd\'^- Hence, p{x) van- 
ishes as x^^^ for a; ^ and as a consequence the above 
integral C exists. 

For fully connected graphs, triangles are the shortest 
orbits compatible with Eq. Q and we have V^^\5 
0) — C, i. e. the probability of scarring does not depend 
on the required quality of the scar. This is in excellent 
agreement with numerical results (Fig. Hfe) and compat- 
ible with the step- like cutoff in ViT) [l9j. In a similar 
way we can describe the statistics of scars on other orbits. 
For example, according to Eq. Q , the square-shaped or- 
bits of Tabled can support scars and we did observe such 
states. However, Eq. (O gives 'P^^\6 ^ 0) 5'^/'^, i. e. 
the probability of strong square scars is much smaller 
than for triangles and consequently they leave no distinct 
trace in ViX). 

In contrast to fully connected graphs, Eq. {7)) allows 
in star graphs scarring on the F-shaped orbits of Ta- 
blelU because the outer vertices have valency 1. Applying 
Eq. ifT^ in this case we get V{6 — ^ 0) ^ (5^^/^, i. e. scar- 
ring is strongly enhanced. As a consequence V{T) is here 
totally dominated by scars, showing a strong maximum 
at X = 0.25 (Fig.[2Jl. It would be very interesting to relate 
this fact to spectral statistics, which for star graphs cor- 
responds to pseudo-integrable instead of chaotic classical 



dynamics |l6l |. 
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